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Improved Overlap Fermions 
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We test exact and approximate Ginsparg- Wilson fermions with respect to their chiral and scaling behavior in 
the 2-flavor Schwinger model. We first consider explicit approximate GW fermions in a short range, then we 
proceed to their chiral correction by means of the "overlap formula" , and finally we discuss a numerically efficient 
perturbative chiral correction. In this way we combine very good chiral and scaling properties with a relatively 
modest computational effort. 



Recent work revealed that the Ginsparg- Wilson 
relation (GWR) § 

{D x , y: 75} = 2 (Dj 5 RD) 

provides the correct chiral behavior of the lattice 
fcrmion characterized by the lattice Dirac opera- 
tor D, if R is a local Dirac scalar [||. However, 
this relation does not imply anything about the 
scaling quality. Our goal is the combination of 
excellent chiral and scaling behavior - as well as 
practical applicability. For details, see Ref. ||. 

In perfect and classically perfect actions, the 
chiral symmetry is represented correctly, since it 
is preserved under block variable renormalization 
group transformations, so it can be traced back to 
the continuum [Q. This is in agreement with the 
fact that such actions solve the GWR. For the 
classically perfect actions, this has been shown 
in Ref. ||, and it is also known that they scale 
excellently. However, they cannot be applied im- 
mediately since they involve an infinite number 
of couplings. Unfortunately this seems to be true 
for any solution of the GWR HQ . A truncation, 
which is applicable to QCD ||, is the "hypercube 
fermion" (HF) with couplings to all sites in a unit 
hypercube. We test the quality of suitable 2d HFs 
regarding: 



• chirality: we focus on the form 2R x ^ y — 
p $x,y (l 1 > 0), where the spectrum <r(D) 
of an exact GW fermion lies on the circle in 
(Fwith radius and center 1/fj, (GW circle). 
We check how well this is approximated. 



• scaling: we test the fermionic and 
"mesonic" dispersion relation. 
Our general ansatz reads D — p v ~f v -\-\, where the 
vector term p„ is odd in v direction and even in 
the other direction, while the scalar term A is en- 
tirely even. In the free case, p v {x — y) {X(x — y)} 
contains 2 {3} different couplings. We consider a 
massless "scaling optimal" HF (SO-HF), which 
is constructed by hand. Its free couplings are 
pi(10) = 0.334, pi(ll) = 0.083; A(00) = 3/2, 
A(10) = 2A(11) = -1/4 (similar to a truncated 
perfect free fermion). The free spectrum approx- 
imates well the GW circle for /i = 1 (standard 
GWR) ||. From the free fermion dispersion re- 
lation shown in Fig. 1, we see that the SO-HF is 
strongly improved over the Wilson fermion. 
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Figure 1. Free fermion dispersion relations. 

We now proceed to the 2-flavor Schwinger 
model, and we attach the free couplings to the 
shortest lattice paths only (where there exist 2 
shortest paths, each one picks up half of the cou- 
pling). Moreover, we add a clover term with co- 
efficient 1, since this turned out to be useful. PI 
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This is not renormalized in the Schwinger model 
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Typical configurations on a 16 x 16 lattice show 
that the deviation from the circle increases with 
the coupling strength, see Fig. 2. 
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Figure 2. Typical SO-HF spectra at strong and 
weak coupling, approximating the GW unit circle. 

As a scaling test we consider the dispersion re- 
lations of the massless and the massive meson- 
type state - which we denote as it and rj - and 
we find again a strong improvement of the SO-HF 
over the Wilson fermion, see Fig. 3. The SO-HF 
reaches the same level as the classically perfect 
action p0| , although it only involves 6 different 
couplings per site (as opposed to 123). 
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Figure 3. The mesonic dispersions at f3 = 6, 
based on 5000 configurations. They are gener- 
ated quenched, but their evaluation does include 
the fermion determinant. 

Hence the SO-HF has a remarkable quality 
with respect to chirality and scaling, but the fig- 
ures also show one unpleasant feature: there is a 



strong additive mass renormalization, which leads 
to uin ~ 0.13 (at [3 = 6). There are various ways 
to move towards the chiral limit: the standard 
method is to start from a negative bare mass, 
but one can also achieve criticality solely by intro- 
ducing fat links with negative staple terms (which 
need to be tuned) || . We now want to discuss yet 
another way, using the overlap formula ]Tl[ |. 

Our ansatz for D obeys — 75.D75. We now 
define A = D — fj,, thus the GWR (with 2R XtV = 
[i 5 x ,y) is equivalent to A^A = fi 2 . If we start 
from some Aq = Dq — fj,, then the GWR does not 
hold in general. However, the operator can be 

"chirally corrected" to A = fiAo/ '\J ' A^Aq. Now A 
does solve the GWR. H. Neuberger |ll| suggested 
to insert the Wilson fermion Dq — D^/ = Ay/ + 1 • 
We denote the fermion characterized by -Djve = 

1 + Aw I \J A w Aw as Neuberger fermion. At least 
in a smooth gauge background it is local (in the 
sense that its couplings decay exponentially) |l^] . 

The overlap-type of solution to the GWR can 
be generalized to a large class by varying Dq 0. 
In particular, if Dq represents a GW fermion al- 
ready (with 2R x y = fi S Xt y), then D = Dq. If 
we now insert a short-ranged, approximate (stan- 
dard) GW fermion like Dso-hf, then it is al- 
tered only little by the overlap formula (with 
fj, = 1), D « Dso-hf- As a consequence, we can 
expect a high degree of locality - i.e. a fast expo- 
nential decay - which is indeed confirmed in Fig. 
4. Also the good scaling behavior of Dso-hf is 
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Figure 4. The locality of the overlap SO-HF vs. 
Neuberger fermion, measured by the maximal cor- 
relation f(r) between two sites at distance r (as 
suggested in Ref. [11]) on a 24 x 24 lattice. 

essentially preserved in the overlap SO-HF, see 
Figs. 1 and 5 for the fcrmionic and mesonic dis- 
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persions. The same is true for the approximate 
rotational invariance ||. At the same time, we do 
have exact chiral properties now (hence = 0). 
If we look at the spectra of certain configurations 
before and after the use of the overlap formula, 
the effect of the latter comes close to a radial pro- 
jection of the eigenvalues onto the GW circle. 
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Figure 5. The mesonic dispersions of different 
overlap fermions at [3 = 6. 

However, in QCD the use of the full improved 
overlap fermion might be tedious due to the 
square root. In view of d = 4 we suggest to eval- 
uate the square root in 

D = n 

just perturbatively around fi. For an approxi- 
mate GW fermion like Aq — Aso-hf this ex- 
pansion converges rapidly, since the operator e = 
A^Aq — /i 2 obeys ||e|| <C 1. It fails to converge, 
however, for the Neuberger fermion. If we per- 
form the perturbative chiral correction, we obtain 
an operator of the form D pxc = fi — AoY. For the 
correction to 0(e n ) the operator Y is given by 
a polynomial in A^Ao/fi 2 of order n. The imple- 
mentation then requires essentially 1 + 2n matrix- 
vector multiplications (the matrix being A resp. 
A'), hence the computational effort increases only 
linearly in n. It turns out that for moderate cou- 
plings already the leading orders are efficient in 
doing most of the chiral projection, see Fig. 6. 
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Figure 6. Histogram of the small real eigenvalues 
of D pxc - based on Dq = Dso-hf - showing that 
the mass renormalization vanishes quickly under 
perturbative chiral correction. 

Therefore, the perturbative chiral correction of 
a good HF combines excellent scaling and chiral- 
ity (and rotational invariance) as well as a high 
degree of locality with a relatively modest compu- 
tational overhead. This method is very promis- 
ing for the extension to 4d, which is currently in 
progress. 
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